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ABSTRACT 


This thesis describes the study and development of a 
workable inverse scattering method for imaging and identifi- 
cation of radar targets. The space-time integral approach 
is used for iterative target shape reconstruction. 

Following an overview of transient electromagnetics, the 
integral equation is applied for thin-wire transient response 
computation. 

The analytical time domain integral equation is derived 
and solved numerically for general conducting bodies of 
revolution. Finally, the algorithm for an inverse scattering 
computer solution is derived and tested under simulation of 


physical environments. 





TABLE OF CONTENTS 


iy INTRODUCTION --------- nner 
BR. OVERVIEW mmm rm 
ERO RACIICAL CONS EDHRATLONS =———=—--<==—------—---- 
Ge SURVEY OF THE LITERATURE ----<=------------- 
bee THESES OBJECTIVES ~-—-—= 99999938 rrr 


II. GENERAL SOLUTION OF THE TIME DOMAIN 
INTEGRAL EQUATIONS -------------------9 7999-7777 


ee THE EMTEGRAM BOUATIONS ===——=>-<<--—--—-------- 
B. NUMERICAL SOLUTION AND CONSIDERATIONS -----~- 


Cs APPLICATION TO THE TIME DOMAIN 
SCATTERING BY A THIN-WIRE -oorr rrr ccc rrr 


fe Le THE TIME DOMAIN RESPONSE OF PERFECTLY 
BONDUGRING BODY <OgMPEVOLUT ION =—— oo pees oo] === 


A. INTRODUCTION rrr mr rere 
B. THE MAGNETIC FIELD INTEGRAL EQUATION -----~-- 
Go) FARVZONE, SCATTERED!) FIEDD =---<-----------"---" 
bee NEAR’ ZONESBAGK—=SCATTBREDSFIELD ==s5=------—— 


IV. NUMERICAL SOLUTION OF THE TIME DOMAIN RESPONSE 
GH THe BODM on REVOLUTION =s==————— oo --—-———-—--s> 


Pen? EDING othe BODYeLNI@seATCHULS (——-—-==s=-——————-— 


B. NUMERICAL SOLUTION OF THE INTEGRAL 
EQUATIONS =-- n-ne rrr 


C. NUMERICAL EXAMPLES ------93 rrr rrr rrr 
V. INVERSE SCATTERING srr rrr rrr 
A. PROBLEM DEFINITION ~----eororr rrr 


ERE ORMULATION OF THE PROBLEM =<=-—----=s=—--—----— 





oF 


D. 


he 


Vie SUMMARY 


oF 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


APPENDIX 


INVERSE SCATTERING SOLUTION PROCEDURE ------ 


NUMERICAL SOLUTION AND EXAMPLES ------------ 


SUMMARY OF RESULTS ------------------------- 


RECOMMENDATIONS ------3- 99-09-2552 -------- 


PenAC EO R RUE RHER ES TUDIESE====——— === = 


A: 


iB: 


E': 


DERIVATION OF THE SCALAR INTEGRAL 
EQUATIONS --------- 


DEMONSTRATION OF THE SINUSOIDAL 
VARIA T LON OH 0 HEm@UiR RE NC SoeWy LU cee ee oeine 


SIMPLIFICATION OF THE SCALAR 
TNTRGRAL FOUATHONSweneee == -- eee 


COMPUTER PROGRAM FOR THE TRANSIENT 
RESPONSES OF ASSTRAIGHT THIN-WIRE =------- 


COMPUTER PROGRAM FOR THE TIME DOMAIN 
RESPONSE AND INVERSE SCATTERING FOR 
THE BODY OF REVOLUTION ------------------ 


RELATIONSHIP BETWEEN RAMP RESPONSE 
AND THE SHAPE ---------------- 9-9-5 


LIST OF REFERENCES ----------- 2-9-2222 = 


eek Or eee OC) TRL ON eT RL i ie ee SS ---------- 


fal 


US 


82 


96 


96 


97 


ay 





LIST OF FIGURES 


Space-time sampling diagram ------------------- 
Geometry of the thin-wire scatterer ----------- 
Current induced at the center of a thin-wire 

by a gaussian incident pulse. (A= 6-10%, 

a/L = 0.0008) --------------------------------- 


Backscattered field from a thin-wire. 
(A = 6-10°%, a/L = 0.0008) --------------------- 


Backscattered field from a thin-wire. 
2 =62.25-107, aft = 090067) =-----—=====------ 


Geometry of the scattering by a body 
of revolution --------------------------------- 


Body contour ---------------------------------- 
Geometry of the far-zone scattered field ------ 
Geometry of the near-zone backscattered field - 
Geometry of the contour sampling -------------- 
Space sampling of whole body ------------------ 
Gaussian shaped pulse incident on a sphere ---- 


Smoothed impulse response of a 1 meter 
radius sphere --------------9 9 nen on -- 


Ramp response of a 1 meter radius sphere ------ 


Measured impulse response of a 30 cm 
diameter sphere ----------- eo 


Geometry of the cone sphere target ------=----- 
Smoothed impulse response of the cone sphere -- 
RalipeeeesDOnSeso: the cone sphere Qe so eee 
Ecomne Cymome Chemcappe d wey bende 


Smoothed impulse response of the 
eappean cylinder -----<=8"------<--------—-——------- 


42 





a. 12 


syeak 


8 


Ramp response of the capped cylinder 


inverse scattering geometry 


Contour estimates of the sphere 
6th estimate of the contour of the sphere 


Contour estimates of the cone-sphere 
to 12 segments) 


(cut 


6th estimate of the contour of the cone- 
sphere (cut to 12 segments) ------------------- 
Estimates of the contour of the cone sphere 
(cut to ll segments) -------------------------- 
6th estimate of the contour of the cone 

sphere (cut to 1l segments) 


Contour estimates of the capped cylinder 


4th estimate of the contour of the 
capped cylinder 


Ramp response of the sphere 


added (S/N = 10 dB) --------------------------- 
Ramp response of the sphere with noise 

added (S/N = 20 dB) --------------------------- 
Ramp response of the sphere with noise 

added (S/N = 30 dB) ----------------- 5-5-5 ---- 
Ramp response of the sphere with noise 

added (S/N = 40 dB) ------------------5-------- 
Estimates of the sphere shape (S/N = 20 dB) --- 
6th estimate of the sphere shape 

(S/N = 20 dB) ------------------- ee ner 
Estimates of the sphere shape (S/N = 30 dB) --- 
8th estimate of the sphere shape 

(S/N = 30 @B) --------------------------------- 
Estimates of the sphere shape (S/N = 40 dB) --- 


80 


81 


83 


84 


86 


87 


8 8 


8 9 


JI 


OEE 


oi 


ye 


215 





io LNTRODUCTION 


A. OVERVIEW 

The classification of radar targets, has become increasing- 
ly important in the past two decades. In military applications, 
the need to distinguish between friendly, enemy, and neutral 
targets, is vital. 

In most radar applications, the only properties that are 
measured are range, angle, velocity and elevation. The methods 
that are used by radars to provide some classification of the 
target generally involve the examination of the echo signal. 
Some of these techniques can be summarized briefly as follows: 

1. High range resolution radar, which uses a very short 
pulse that can provide sufficient range resolution to obtain 
peeOugHeprotilewof*™the illuminated part of the target shape, 

2. Cross section fluctuations analysis which can provide 
information on the size of the target. 

3. Synthetic aperture radar that gives a high angular 
resolution. This method is typically employed in ground 
mapping. 

4. Analysis of the backscattered polarization from a target 
can provide a means for discrimination between types of 
lea Gets. 

5. Use of nonlinear effects due to junction of metals. 


The backscattered spectrum is analyzed and checked for 





harmonics appearing in the received signal, that are typical 
for certain type of targets. 

6. Doppler modulation by propellers, the rotary wing in 
helicopters or even hull vibrations in armored vehicles. 

None of the methods cited can be used to obtain all the 
information on the target shape and size, they only can be 
used to discriminate to a certain extent between different 
types of targets. Other methods currently used involve 
passive reception of radiation by targets. 

Elint receivers or IR detectors can be used to monitor 
radiation from targets, and, based on previous intelligence, 
or Signature catalog, associate a type of radiation to a type 
Smee Garget. 

The inverse scattering technique,that is the subject of 
study in this thesis, uses the basic properties of the transi- 
ent response to an impulse excitation (impulse response) in 
order to obtain information on the size, shape and also material 
composition of the target. 

In principle, the smoothed impulse response can be ob- 
tained by measuring the backscattered fields at all frequen- 
cies up to some upper limit. In most cases this would be 
impractical so, instead, the target is illuminated by a very 


short pulse (with a very wide band frequency spectrum ranging 
from the low frequencies in the Rayleigh region Up te higher 


resonance region frequencies). 
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The impulse response contains all the information needed 
to characterize the target. The higher frequencies give the 
information about the structure. This is related to the 
optical theory of diffraction which is a good approximation, 
under certain conditions, for shape prediction up to the 
shadow boundary. Physcial optics, however, fails to predict 
the phenomena appearing in the shadow region, isuch as*the 
creeping wave which is due to lower frequency waves travelling 
around the body. 

A different approach that has been used to retrieve tar- 
get configuration information is related to the singularity 
expansion method (S.E.M.). This method seeks to classify 
scatterers through the complex residues and poles of their time 
domain signatures. Briefly, in this method the transient E.M. 


response of a conducting body is characterized as a series of 


complex exponentials (or damped sinusoids): [i], [2] 
N sat 
ce) oe 
m=1 


where Rn is the amplitude (residue) of each mode of complex 
frequency s_- The true function has a complex frequency 


representation of the form: 


N R 
iS) ace 


m=1 


MT 
SI > 
mM 


The poles Sa depend only on the geometry (independently 


of the orientation), and the residue Ru depends on both 


Ii 





excitation and geometry [3]. Thus, once the poles and residues 
are known, they provide a means for target classification. 
The problem is to find the poles and residues. 

A method is suggested using Prony's algorithm in [2] for 
extraction of poles and residues from time-domain signatures. 
The method has been tried successfully on wire geometries and 
Simple shapes. 

Another important application of the impulse response of 
a target is the study of its radar cross section (RCS), since 
the Fourier transform of the impulse response, (the frequency 
response) is directly related to the RCS by the following 


relationship: 


where 


and 


PL (t) is the impulse response of the target. 


Finally the time-domain analysis of the impulse response of 
a body is of great importance in the study of the electro- 


magnetic pulse (E.M.P.) response of targets. 
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B. PRACTICAL CONSIDERATIONS 

The practical study and measurement of the impulse re- 
sponse involves the building of a time-domain scattering range. 
The availability of generators having very short pulse-widths, 
fast rise time, and high voltage outputs, together with fast 
sampling oscilloscopes and small computers makes the task 
feasible. 

For the solution of the inverse scattering problem the 
time domain scattering range developed at the Naval Post- 
graduate School [4], was employed to measure the impulse 
response and demonstrate the validity of the theories developed 
in this thesis. The transient response data measured by using 
pulse generators with very fast rise times can validate the 
calculation of the time domain integral equation solution. 
More important, it provides a vehicle for testing various tar- 
get imagery and I.D. schemes, as well as investigating broad- 


band RCS reduction. 


See oURVEY OF THE LITERATURE 

One of the first efforts to determine physical properties 
of electromagnetic scatterers from time-domain analysis was 
due to Kennaugh and Cosgriff [5], who employed the physical 
optics approximation to calculate the approximate backscattered 
impulse response of a flat spheroid. 

A summary of this work was reported by Kennaugh and Moffatt 
[6], who extended the physical optics approximation to transient 


response calculations. At this stage the computation of the 
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impulse response was approximated by using the physical optics 
approach or the inverse transform of the frequency domain solu- 
tion. Then, in 1968, Bennett derived in his Ph.D. thesis [7] 

an approach for solving the time domain integral equation 

Hates cab by using a time-stepping method. Bennett used a 

form of the equation called the magnetic field integral equation 
(or M.F.I.E.) for two and three-dimensional surfaces. 

Sayre and Harrington [8] used the same time-stepping tech- 
nique to derive the transient response of a wire, but they 
employed another type of equation called the electric field 
integral equation (EFIE). 

Additional work in time-domain electromagnetics has been 
reported in several more recent papers as, for example, the 
report by Miller, Poggio and Burke [9] on time-domain analysis 
of thin-wires. Summaries and review of the time-domain electro- 
magnetics are due to Mittra [2] and [10] and Miller [ll]. 

A comprehensive review on the subject and its application 
has been published by Bennett [12], which gives the reader 
an introduction and initiation to time-domain electromagnetics. 
Literature dealing with the inverse scattering problem is due 
to Young [13] and Shubert, Young and Moffatt [14], which des- 
cribes the image generation of a target from the ramp response 
waveform. 

Moffatt and Mains [15] describe a method of detection and 
discrimination of radar targets using the idea of ramp response 


but with multiple harmonic radar frequencies. Finally,a very 
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recent report by Bennett [16] gives the general approach to 

the time domain inverse scattering solution. The method sug- 
gested in that report is being used in this thesis to approach 
the problem of transient target imaging. This list of efforts 
is by no means exhaustive, as many other authors have published 
papers dealing with transient electromagnetics and target 


imaging. 


D. THESIS OBJECTIVES 

ahe objective of this effort is to develop the analytical 
and numerical solution of time-domain inverse scattering 
problem. It follows the development of the time-domain 
scattering range at the NPS [4]. 

The main effort has been conducted toward the formulation 
and numerical solution of the time-domain integral equations 
for conducting bodies of revolution. 

A computer program is developed to solve for direct and 
inverse scattering. The algorithm was used for experimental 


display of the shape of the target under test. 
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Il. GENERAL SOLUTION OF THE TIME 
DOMAIN INTEGRAL EQUATIONS 
This chapter is a review of the development of time domain 
integral equations, their solutions and the associated numeri- 
cal considerations. Detailed development of the equations 


Sammbeatound in {2,10,ll1]. 


A. THE INTEGRAL EQUATIONS (E.F.I.E. and M.F.1.E.) 
The basis for derivation of the equations are the general 


time dependent forms of Maxwell's equations: 


> > 0 > > 
Ween G,c) =. = =e Ho pla aie. 
FF @,s) = = ce, E(Z,t) + F(Z,t) 
Ve he E (r,t) = p(r,t)/e, 
_—> 
Ve Boece) = - 


Starting from these equations, and using the expressions for 
the magnetic vector potential, A, and the electric scalar 


potential, ¢, which are related to the fields by 


> > 
5 (G9 45, 


lI 
] 
J 
<> 
ry 
ct 
ay 
ct 


E(r,t) 
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one can obtain the expression of the magnetic field integral 


som@icrom (Merl; E.) for a perfect conductor [2] 


Jo(r,t) = 2nxH (ie) 
ibaa i) pes ge R 
+ ae ned le a7 t+ RY, (eT) “ee dos (2 els) 


e = u 
eC > a n Od ae 
nxE (r,t) = se 8 lage ree OM ao) 
S 
(= + = 57) R]ds' (2.2) 
EGR 


In equations (2.1) and (2.2), 


n is the unit vector normal to ENe@wCcOnductor wae 
the observation point r. 


me * * e 
J_(r',t) 1s the current density at the source point r! 
at a retarded time T, 


o(r',t) is the charge density at source point atthe 
retarded time T, 


Here t iS given by 


| = t- R/c 
and 


. + Te 
R is the distance from the obServation point r to 
ShemsouEce point Giyme.cuw 


-- > ee 
R = re-fr 


ey 





ze1gyS eee = TGS : mae 
E a,t) and H (r,t) are the incident electric and magnetic 


Fields at the observation point. Because of numerical con- 
straints, the equation most suitable for closed surfaces is 
the MFTE equation (2.1) while for thin wires and open struc- 
tures the EFIE equation (2.2) is used. 

Equation (2.1) states that the current induced on a point 
r of a body at a time t is due to the incident field Rone at 
that point and the sum of scattered fields there from earlier, 
more distant locations of the body. The interactions between 
the samples of the body are displaced in time by an amount 
equal to the transit time of the fields between them. 

In the integral, the observation point (R = 0) is excluded 
Since this "self-term" is included in the incident field por- 
tion of the equation. This determines the general solution 
for the MFIE equation which is solved directly by marching 
on in time. The solution is started at a certain time (t = 0), 
and due to causality the integral part of the equation will 
be zero Since there are no past currents. Then, as t increases, 
the integral is solved with previous currents already computed. 
This technique is a time-stepping procedure for solving initial 
Yalue’ problems. 

Kor Sthes brah ,@equation (2.2)> the§method ais?slightly differ- 
ent but uses the same time stepping procedure. The solution 
Sethi smecuatiton ws applied to the problem of E€ransiene 
Beatteerimo byeamthin wired, #ln the next chapter. 

ijeesomuciOn of equations (2.1) and (2.2) gives the cur- 


rents induced on the body by the incident field. Once the 
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currents are known, it 1S a straightforward step to compute 


the scattered fields. For the far-field case the field is 


elvene py : 
eee ty E oe ate ty xe aS" (2.3) 
o! 4mroc F OT aes “R ; 


where an is the unit vector from the source r' to the far 
observation point and ty is the separation distance. ay Galo) 
is the current at the source point at a previous (retarded) 


time T = t-xri/c. 


B. NUMERICAL SOLUTION AND CONSIDERATIONS 

The solution of the integral equations implies finding 
the currents induced on the body once illuminated by an inci- 
dent field. Since the solution is to be carried out numerically, 
the first step is to represent the currents by their space and 
time samples Ig 9 Nini) in a discrete manner. Ty 2 Weiae is the 


Eh h time 


current at the 1 segment AS (or patch) at the ;° 
increment AT. The total current will be approximated in terms 
of all samples at all times and spatial segments via the basis 
function expansion: 
No Nw 
cet) = 2 ay J; 4 (s/t) U(t,) V(s,) (2.4) 


where No is the number of space samples on the body while Nin 
is the total number of time increments, and the pulse basis 


functions are: 


Joe, 





tl at the center of the jth time interval 


S22) a 

0 elsewhere 

1 at the center of the ith space patch 
V(si) = 

0 elsewhere. 
The space time sampling is shown in figure 2.1. The 


sampling of the currents already suggests that the body is 
divided into segments (or patches) AS and the time is sampled 
at AT intervals, while the total observation time goes up to 
Nx x AT. 

It 1S very important at this point, before attempting to 
solve this equation, to define the basic criteria for proper 
selection of the space and time sampling as well as their 
relationship to the shape of the exciting pulse and the body 
size. 

The incident field is chosen here to be represented by a 
gausSian impulse due to its rapidly decreasing features and 
also because most of the short pulse generators will produce 


such a waveform. 


The gaussian impulse is given by 


~—A t 
g{t) = e 
Its amplitude falls to 1/10 of its maximum value at 


oe = 1.5/A. Its frequency spectrum is given by [2] 


PLP 2 
mia r(git)} ote" = | 
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computation point 
space |j 
j+3 
}+2 


j+l1 





integration region a 


Figure 2.1. Space-time sampling diagram 
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This response reduces to 1/10 of its maximum value at 
fax ~ 0.5A. The time sampling AT must exceed the Shannon 
rate to avoid aliasing problems if FFT's are employed later, 


1.e., 








The spatial sampling AS must be chosen such that 
AS > ATxC (2.6) 


otherwise there will be coupling between the currents of 
adjacent patches. Violation of (2.6) will also mean that the 
stepping procedure for solving the equations, which assumes 
that during one time interval the currents at one patch are 
not influenced by adjacent currents, will no longer be satis- 
fied, thus producing an incorrect solution. 

Another consideration is that the space sample points must 
be close enough to resolve adequately the spatial variation 
mrcthe inerdent field over the body. From Eq. (2.5) and (2.6) 


we have: 
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oe 
As > —mi0 


This latter relationship imposes some minimum ratio of the 

size L of the body to the minimum wavelength (or equivalently 
the minimum pulse width). The ratio L/r nin determines the 
resolution and the usefulness of the time domain analysis, 

and it must be kept high enough to maintain adequate resolution. 
If the body is to be represented with good fidelity with a 


certain minimum number of samples No we will need 


Ng an 


Le 5 


This means that if dpe increases (1.e., wider pulse) L must 
increase or equivalently for a small size target Mena must 
decrease (1.e., Shorter pulse). 

Having presented the equations, their general solution and 
the requirements for the various parameters involved we can 


move on to the implementation. 


C. APPLICATION TO THE TIME-DOMAIN SCATTERING BY A THIN-WIRE 
The first application that was studied was the computation 
of the impulse response of a straight thin wire. This example 
was chosen since it gives a good understanding of the physical 
processes and the numerical techniques used in the solution 
of time domain integral equations. The basis for this compu- 
tation is given by Sayre and Harrington in [8]. 
The problem is formulated as follows: a thin-wire of 


length L and radius a aligned with the z axis is illuminated 


pe 





by an impulse field of the form 


BINS (4) = i max 
The basic formulas are in terms of the retarded potentials: 
The boundary condition : 

inc 


E x: 
Saye Se BS: 


where the retarded magnetic potential is given by: 


t 
Aj(z) = gf ee ae 
Ax1s 


and the retarded electric potential is 





Pr il q(z',T) 
o & aire J R dz 
@ Axis 


where the continuity equation between charge and current is 


CEU, Sr 


Z, ot 


In these equations I(z',T) and q(z‘',t) are the current and 
charge at source point z' at time T = t - R/c, while R is 
the distance between the source and observation points. The 
geometry 1S described in figure 2.2. 

Using the numerical formulation expressed in [8], and 
applying the time stepping procedure, a computer program was 
written (Appendix D) to compute the scattered field from a 
thin wire of length 1 meter and radius 0.8 mm. The wire was 


divided into 40 segments and the sampling time AT = AZ/c = L/40c. 
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Figure 2.2. Geometry of the thin-wire scatterer 
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The incident field was taken to be similar to the pulse 
transmitted by the generator of the NPS transient scattering 
range, with linear polarization parallel to the wire and 


incidence normally on it. 


. ~a* (t-t___)* 
Lc max 
E = ea 
A = So ip gece 
2 = 18-AT (about half transit time along 


the wire) 


The program computes the currents at each spatial segment with 


up to 400 time samples. This corresponds to 10 transit times 


across the wire. The current at the center of the wire is 
given in figure 2.3. The backscattered field is given in 
figure 2.4. 


No significant change was obtained for 20 segments and 
200 time increments. Comparison of this computation with lab 
measurement was done by Hammond [4] and very good agreement 
was found. 

Similarly, comparison of the computed results found in 
[ll] with a ratio of L/a = 0.00667 and A = SB Ie gave 
exact agreement. The backscattered field for this case is 


shown in figure 2.5. 
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Figure 2.4. Backscattered field from a thin-wire 
(A = 6°10”, a/L = 0.0008) 
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til. GIMME DOMAIN RESPONSE OF PERFECTLY 
CONDUCTING BODY OF REVOLUTION 
A. INTRODUCTION 

This section considers the solution of the time domain 
response for a body of revolution, when the exciting wave is 
propagated along the axis of symmetry of the body. 

The use of the axial symmetry of the body will result in 
considerable simplification of the computation of both the 
induced currents on the body surface and the field scattered 
from the body. The algorithm developed in this section will 
be used later for the solution of the inverse scattering 


problem. 


B. THE MAGNETIC FIELD INTEGRAL EQUATION 
In [7] Bennett derives the integro-differential equation 


for three dimensional scatterers: 


> “A 1L.> 
UtGeejee — wax HH (r,t) 
i, i ree ln oo 
me en Se ee eae GD) 
A R 

where: 

T = +t - R/C (retarded time) 

r is the observation point 


KY 


is the source (or integration) point 
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Ree rea 


A is the surface of the scatterer 
aS is the incident magnetic field at the observer's 
pointe 
The geometry for the body of revolution is given in figure 
ol. 

For the actual problem the incident magnetic field is 
propagating along the z-axis (axis of rotation of the body), 
and is directed towards the positive x-axis. the body lies 
completely to the right of the positive z-axis, and is defined 
by the continuous contour function, p(z), which is the radius 
of the body at a point 2z. 

The prime coordinates denote the source or integration 
point, while the unprimed denote the observation point for 


which the current iS desired. At each point on the body, there 


) 


wN ~s 


is defined the triad of orthogonal unit vectors, (a_.a ya 


s' 9 
where: 
a: unit vector normal to the surface 
ay? unit vector tangent to the surface and directed 
towards the positive ¢$ direction 
am: unit vector tangent to the surface and directed 


along the curve representing the longitudinal 
contour of the body. 


For the following derivation we will define the angle a (or 
a') as that sustained by a. (or a‘) Withothemz axis, co (or a) 


being poSitive when a. (on a!) points outward from the Zz axis 


SJL 





and negative if it points towards it. The angle $ (or 9') 
is the angle defined by op (or o0') with the (x,y) plane. 
Eguation (3.1) states that the current density at the obser- 
vation point r at time t is due to: 

1. The current induced by the incident field at (r,t). 
This part represents the physical optics approximation, 


>i o> 
2a, ate (te) a 


2. The current due to the flow of currents at all other 
points of the body at a retardated time Tt = t - R/C, where 
Ris the distance from the integration point to the oberva- 
—i10n point. 


According to the geometry described in figure 3.1 we can write 


2 


R = (0 +p’ spp" com(o' <6) + (z" -z)*)1/4 


The current on the body will flow only in the S and $ direc- 
tions and can be expressed (for both primed and unprimed 


coordinates), as 


=e: A A 
Deby ft) a J(0,2,t,9) ae ear), (Oy 27 © Ad) a 


p p 


where J (p,2,t,$) and J,(9,z,t,%) are the scalar values of 


p 


the current components and are functions of the position on 
the scatterer. We will denote them simply by J and J SO 
that 


BZ 





Leu e, sie 





Geometry of the scattering by 
a body of revolution 
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Eq. (3.1) can be written in terms of its components in the 


s and ¢ direction, 1.e., 


J,(rrt) = (2a, x (r,t) h, ‘aq Jazlgt +e Fella, x5" «Rl aa’ 
3.2) 

“ ey 1 lee eo eee , 

J,(r,t) = [2a, xH (rt) - an Jak tS 57) lay 9 sla hsp 
G3) 


In order to obtain the scalar expressions of the two last 
equations we have to perform the vector multiplications: 
This is done in Appendix A and as a result, equations (3.2) 


and (3.3) become: 


y aa et Eig A wre: , 
i = Ht (z,t)sin $d +5— | Sale tos | [OR J 4V sin 8]dA 
AR 
(3.4) 
—_- 1 fii?  l 
- le 2H* (2, t)cos ¢ sin a+>— waele py TRI [IZ U sin B +d Fl IGA 
(3.5) 
where: 
oe OF 
Py = (z' -z) cos8 Sind + (9p -—o' cos8)cos a' 
V = z2z' -2 
F, = (op cosg - o')cosa + Sinacos 8(z' -2Z) 
U = (z' -z)Sina-:sin a' +o cosa sina' -p'sina-cos a' 
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We note that the incident field H~ is dependent only on z and 


t (independnet of ¢ due to symmetry). We also note that: 
aA =8 —* do' ds’ = p* aB ds? 


The integrals in Eq. (3.4) and (3.5) can then be divided into 
two parts: integration over 8 from 0 to 2m, and integration 
over ds' from 0 to S (S being the contour curve as shown in 


megure 3.2). 
p (Zz) 


L Z 
Figure 3.2. Body contour 
The next step is to obtain further simplication of the inte- 
gral equations (3.4) and (3.5). 
The currents in these equations are functions of their 


position on the body: 


Cy 
ti 


J5(2,0,6,t) 


Cy 
t] 


6 ee 


But we remark that the currents due to the incident field have 


a Sinusoidal variation with ¢, as follows: 


J = H™ (z,t) cos ¢ + sina 


p 
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J. = H- (z,t) Sin 9¢ 


so that we can assume that the currents will have the same 


respective relationships with $6. 


J5(Z,0,t,$) (2 eee isin On J.:sin d (3.6) 


S 


J,-cos 6 O37) 


UZ, 07 acoso j 


Jy 6210 1t,9) b 


This assumption is demonstrated in Appendix B, using Maxwell's 
equation. 

After substitution of equations (3.6) and (3.7) into the 
mecegrals (3.4) and (3.5), the integrals over the angle § can 
be simplified as shown in Appendix C, so that the scalar 


integral equations reduce to their final expressions: 


S 27 


es T 1 , 1 dk iE | 
2 = 2H (z,t) + a f o'ds f acme sal 
0 0 R 
. ees 
-(JLF, cos 8 +Vsin oie Mere (3 36) 
. S 27 
Ue = 2H (z,t)sin at T= f 0 ds f alate eae 
0 0 R 
‘? 
- 1 ' 
[Je U sin sg + F. cos 8 glee eo 9) 


in which the currents Jax and J (Jo and J!) are only functions 


? p 
of 0, z, and t (9',z',t) and independent of the angle 94. 
The significance here is that the surface distribution of 


the unknown currents is reduced to a one-dimensional variation 


of current along the generating arc of the surface of revolution. 
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This will allow considerable reduction in computer time and 
storage when carrying out the numerical solution of the 
equations. 

Once the currents along the contour of the body are known, 
the currents at any other point (at some time t and some z 


position) will be given by 
Jg(0,2,t,6) — J,(0,2,t)-sin 0) 


J, (0 ,Z,t) “cos d. 


Diy OA TR 


C. FAR-ZONE SCATTERED FIELD 

In most cases we are interested in the far scattered 
field. At a large distance from the body the scattered field 
1s given by: 


1 


s.> 0 > > A : 
H (rt) Fire Jap (Flt) x ap lda (C370) 


where 


T is the retardated time t - R'/C. 


R' is the path difference between the integration 
point to the observer and the origin to the 
observer. 


S(r',t) is the current at the surface integration point 
retarded in time, while a. is the unit vector 
pointing from the source point towards the obser- 
vation point. The geometry is described in 
@igure®3.3. 


r is the distance from the origin to the observation 
point. 


oy) 









direction of 
propagation 


a 


A is the projection of the source point on the 
plane (x,2Z) 


Figure 3.3. Geometry for the far-zone 
scattered field 


38 





In the following, the scattered field will be derived 
for an arbitrary elevation angle §@. Due to symmetry with 
@ there is no loss of generality if we observe the scattered 
field in the plane (x,z) where o> = QO. 


According to the geometry: 
Reem —iz' cos, 8 + 6° cosid' sin 6] 
foe tt — R/C; 


R' is negative whenever the observation point is closer to the 


source than the origin, and positive otherwise. 


a = Sin 6 a + cos 8a 
R x Fi 


Performing the vector multiplication in equation (3.8) and 
separating the field into its cartesian components, gives the 


followoing set of three scalar equations: 


= =, 1 t 2 t 
Hy, (8,t) = imc G Ji a7 lS Sin 24 Sin a' ae © d'|Jcos 6dA 
(35.9) 
H°(68,t) = : J3 a l(J, -d_’sina')sin $' cos ¢' cos 6 
ad 4mr_ Cc of oT 0) S 
+ J. sin ¢' siné cosa']dA' (32 10) 
Ss rs v = ' t 
ac =) =) 4 ¢ . a | sin @(J_sin* d'sin a’ +3 ,cos Say ]da 
(See) 
After substitution of dA' = o'do'ds', and using the odd 


properties of some of the terms appearing in the integral on 


oo 





o' (as it has been done previously in Appendix C), the three 
scalar equations representing the components of the scattered 


field reduce to the following expressions: 


HY (6,t) = cos 6 H® (t) eres 
Hy (6, t) = - sin 6 He (t) (3.13) 
aay 2 @ (3.14) 
ane } 


where the common term He (t) 1s given by the following integral: 


1 S 27 


S ad C) me 2 
HA (t) = tz oc = o'ds' y) a7 (I, Sin ?' sin a’ 


+ J, cos*o') dé! (3.15) 


Equations (3.12) to (3.15) represent the scattered field for 
an arbitrary elevation angle 68. The special case of interest 
is the backscattered field where §6 = 7. For this case there 


is only one component of the scattered field, which is 
S = TS 
H(t) = H(t) 
and the retarded time 1s: 
tT @=8 E =<R'/C 


where 
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All the terms appearing inside the integral in equation 
(3.15) are independent of ¢' so that the integration on ¢' 
can be carried out analytically. Doing this yields the 
following expression for the backscattered field: 


S 


il. 3 ; 
J 57 bd (e',2',T) sin a' +7 


S = me ' ' ' ’ 
H(t) = az c (O32, tee ds 


> 
er kG) 


In equation (3.16), the currents J. and J, are those given 


p 
by equations (3.8) and (3.9) respectively, used at time 
T=t- 2/C (2 > 0), where z/C is the retardation time rela- 
tive to the origin. The backscattered field equation is 
independent of ¢ which means that we can integrate the time 


derivative of the currents, along a Jine on the contour of the 


body. 


D. NEAR ZONE BACKSCATTERED FIELD 

In the Time Domain Scattering Range developed at NPS, 
the distance from the antenna to the target under test is not 
long enough to be considered for all frequencies of interest 
as a far region. Consequently, it is needed to develop a 
formulation of the scattered field in the near zone. For this 
case the geometry is described in figure 3.4. 


The general equation for the scattered field is given by: 


Qo 


yt) = ge flptge seer) «RIS (3.17) 


AOR OT 
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direction of 
propagation 
”) 


observation 
point 


Figure 3.4. Geometry of the near-zone 
backscattered field 
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where 


4 
il 


t - R/C, and 


qaA' 


o' do' ds! 


12 12 


2 2 
R¢“= (ro + z'") + ¢ 


Peter carrying out the vectorial multiplication, dividing the 
field vector into its cartesian components, again dropping 
all the odd terms in $' and integrating analytically on $', 


equation (3.17) reduces to: 


. 


ne 


ss 


S = | ee 
pet rt) Fae J nT +e 


P{d_(e! cosa't - (ro +2z')sin qa! 
0 R 


Lom, 


1 


eer + 2) yO OS 
ee, (3.18) 


Equation (3.18) is the expression of the near backscattered 


field where: 
Rese (ir +2") + 6°7) 


The currents J(0',2',T) and UP re er are given by equations 


(3.8) and (3.9), being used at the retard time tT. 


T = t- (R'=-r)/C (referred to the origin) 
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IV. NUMERICAL SOLUTION OF THE TIME DOMAIN 
RESPONSE OF THE BODY OF REVOLUTION 
In the previous section the analytical equations for the 
currents and scattered field were derived. The equation ob- 
tained will have to be solved numerically, and the first step 
will be to divide the body into patches of approximately equal 
curvilinear surface area. This is followed by the setup of 


the algorithm for the solution. 


A. DIVIDING THE BODY INTO PATCHES 

The geometry of the body is completely specified by the 
memction 9(Z), (0 < z < L); which represents the body contour 
function S. Since the integration is to be carried over S, 
the contour will be divided into NS segments of equal length 
AS. Consequently, the body contour function is fully defined 


by NS+l points. 
P, (2) 2 = 1 to NS+1 


The space samples of each segment AS is at the middle of the 


segment and is given by 


Pp, = (Po, + Py) /f2 

g =e = 71,..2.,NS G4, 1) 
ee Syme Ao) 7 
The angle a for each sample point is the slope of the curve 


at that point and is approximated linearly by 
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p ae 
Ose = Arctan (th) - Ca Lees, NS ae) 
¥ Q+1 g 
Figure 4.1 shows the geometry of the contour sampling. Each 
segment AS defines a ring at Zs with radius Pas All the rings 


will be divided into patches of approximately equal surface 


area, GA. , according to the following formulas; 


= ® = 2 
da, = DS = pp. seABe ao AS 


where AB. is the angle between two adjacent sample points on 


the same ring. 


This suggests a division of each ring into an even number of 


patches NPI which will be 


iO. 
NPI = 2 x INT (==) (4.3) 


The function INT means the integer value of its argument so 


that the increment angle for the ae ring becomes 
" 2 Ti 
48; > NPT ae 


Thus, following this procedure, any arbitrarily shaped body 
of revolution can be subdivided into patches of approximately 
equal surface area. Figure 4.2 shows the geometry of the 


Space sampling on the whole body. 
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Figure 4.1. Geometry of the contour sampling 
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‘oe space sample point 


Figure 4.2. Space sampling of the whole body 
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B. NUMERICAL SOLUTION OF THE INTEGRAL EQUATIONS 
The equations to be solved are given by Eq. (3.8) and 
(3.9) in previous sections. 
For the numerical solution the currents will be represented 


by their space-time samples: 
Mo, 2, c) = Jute, 


where i is the ie space sample on the contour function and 


nis the ee time sample 


a =) tee. .,NS 
Hime otro. NL 
t = n-AT 


The time being sampled at a rate AT. Similarly the incident 
field is sampled and will be denoted by its mth time sample 


at the +th segment: 


wi (z,t) = H*(i,n) 
The general procedure for the numerical procedure was already 
mentioned in Chapter II, and will be reviewed in more detail 
here. The incident field, H, is assumed to be zero at all 
times less than some reference value (which will be set equal 
tolAt). Hence the total current on the surface of the scatterer 
is also zero prior to time t° 

At the next time interval t, = tot At the incident field 


iL 


reaches the scatterer at z = Zy (first space sample), thus at 


48 





time ty the current will be simply twice the tangential com- 
ponent of the incident field. As time marches on, the incident 
field covers more of the body and the current will be due to 
the incident field plus the fields due to currents at other 
points on the scatterer at earlier times. So, by marching 
on in time, the currents at each space sample, for each time 
increment, can be computed. 

The currents appearing inside the integrals in equations 


(3.8) and (3.9) are the retarded currents, and are given by 


Ya Pay he ot = Pe 8 08 eS, 


The point (o9',z') is represented by the k-th space sample 


while the time Tt will be: 
T = t - R/C = neAt - R/C = gQAT. 


In the general case tT will not be an integer multiple of 
AT, so that the current cannot be represented by any of the 
previously computed samples. In order to find the value of 
the current at time tT, there is a need to interpolate between 
the sampled known, discrete values. 

For a better accuracy use was made of the four point 
Lagrangian interpolation formula using two samples at each 
Ct 


Side of the value gAT if t t, are the four samples 


Te ea 


used for interpolation, such that: 
J t. + At = t, + Zh t, + 3At 


te = gy, Me 
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3 
to = 99> At 
Tl = 9) its 


Gi" Gor Gr Gy are integer values, and 
c= age At (g real number) 

such that 
t < te < T < t. = ty 


The interpolation formula using those four points will be 


given by: [17] 


4 f i= ~t; 
Jeog(t) = PARA 2-9 '*n)! (4.5) 
if¥n 


This interpolation formula must be adjusted in two cases: 


1. Whenever t, or t, 1S greater than the actual time 


3 4 
of computation t (in order not to interpolate in 
the future). 

2. Whenever ty or ty are less than the beginning time 


of computation te (due to causality). 
The time derivative of the currents will be simply expressed 
as the derivative with respect to t of the interpolated 


value, which will give the following formula: 
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UP ge ’ 
ome 
Defining: 
Ti = g 
qT. = g 
qT, = g 
Ty = g 


The interpolation 


derivative can be 


oes he a T 


Js g'tr9) : 





; 5 4 tot, 
arl I (;-)J a (4.6) 
n=1 - i=l La eh Sty is 
1) 


formulas for the current and its time 


written in their numerical form as follows: 


oes oo J 


2 2 5 4 


se Fs,60 194) - 6 


eT 


an, a ——,-—~ Jee.) 


"4 


and 


DJ = 


S,0 OT 


<— eo Ree 


is wget td) TT ae Vek 


Sh ae diel 


12 (4.7) 


FE 


2 2 1 ee. 


5 
6AT 


Sey tte S 4) 


eee Lek a 
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3° 4 2 4 


ia Git mo a, Ge oy! 


le tT ae a 


, 21 


2AT 


A°l 3°74 


Me pag Ma we (dW Mi cid ee 


4 1 


2AT 


2 2 4 
Je SORey (4.8) 
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Now we can write the numerical expression for the currents 
in equations (3.8) and (3.9) by representing the continuous 
integration by numerical summation, and expressing the cur- 
rents in terms of their space-time samples. 

For simplicity, we will write first the expression for the 


integral part of equation (3.8) and (3.9) which will be 





denoted by Jas! and Too: 
, : 7 
J lc oid) = Aa AB OD dL 
Ss Pa8 = 7 ae 1 2 m= 52 


m41 for I=k Rim 


VesI* (m) Jy (irg) +F 
-[ 
“Seta 
W4sT—(m) DJ (i,g)+F,CO(m)Dd .(1,g) 


® 
i = ] (4.9) 


1CO(m) Tig) 


NP. 

= 1 i 

moe) 5. oF 2 AB TP ASy 2 
i=l m=L co 

m#1 for i=k 





U-SI* (m)J, (ig) +F,CO(m) J, (ig) 
- [8 + —___2 
kim 


n = (4.10) 


where: 


J g (4-9) are given by equation (4.7) and DI. aes, 
4 by equation (4.8). ‘ 


Z 2 2 
Ro = (of+p% - 20, 0,CO(m) + (z, -2,) yo 


kim 


CO(m) = Cos{(m-1) AB, ] 
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SI(m) = sin{(m-1) 48, ] 
a — a; = Z 
U = V sin Oy Sin Os + Pp}, COS a, Sin a, 72; Sina x COSA 5 
Fi = V-CO(m)-sina, + (p,-9;,CO(m)) cos Os 
FE. = cos by (9, cos (m) - p;) +V sin a. CO (m) 


NS represents the total number of segments 


NP represent the total number of patches at the 
TEA yimg. 


The total currents can now be expressed in terms of Jas and 


Tog! which yields: 


cf Mice 
J, (k,n) = 2H ee + Jog (Kn) (4.11) 


J, (k,n) Os Bia) (4.12) 


> 


a , 
-2H (2, ,n) x Sin (ay) +S a6 


J. and J, are the values of the current components stored and 


p 


used for computation of the backscattered far field, which 


will be expressed in numerical form as follows: 


NS 
S ee! ; ; ; 
Mm (n) = tz .c =, [DJ _(i,g) sina, +DT, (1,gIAS 505 (oaks) 
where 
yA 
ne 
J Cite 
and DJ. g ft°9) are given by equation (4.8). 
Sf 
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Equations (4.7) and (4.8) used for interpolation of the 
currents and their time derivative are given for the general 
case where: 

i” or > 1 (after the incident field had reached the 

specific point i on the scatterer) 

Di G, <n (so that no interpolation in the future occurs). 
In cases where these two conditions are not met the interpo- 
lation must be adjusted to include only those points which 
satisfy causality with no interpolation in the future. 

The numerical formulations derived in this section were 
programmed in FORTRAN to compute the currents and the back- 
scattered field, due to the excitation of a body of revolution 
by an impulse type field (or a ramp). The input data to the 
program is the function 0(z) along with the length of the 
Spatial segment AS and the time increment AT. The program 


1s given in Appendix E. 


C. NUMERICAL EXAMPLES 

The first example to be considered is a sphere of radius 
1 meter. This special example was chosen since its result 
appears in various publications and comparison can hence be 
made. 

For the present example, use is made of the guidelines set 
in Section (1.B) for the choice of the incident field, the 
Sampling rate and the space sampling. The incident field for 
the case of the impulse response is represented by the gaussian 


shape 
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a 2 
=e) (b=t -zZ/c) 
4 max 2a AG 
@: 
H*(z,t) = 
0 ; otherwise 
where: 
A = Roma age 
which yields 
AT = Q0.2/C sec 
mac pulse is described in figure 4.3. 
X 
<—— 2 a/c : idee Ne 


Figure 4.3. Gaussian shaped pulse 
incident on a sphere 


The pulse reaches its maximum value at the center of the sphere; 


(ai = a/C. 
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The spatial segment, with 15 segments 


AS =e = FO. 209 m < ~ CcCAT 


The number of time increments observed was such that it 
covers 6 times the transit time across the whole sphere (60-AT). 
The results obtained for the backscattered field are given in 
figure 4.4, on which the results of Bennett [7] are superimposed 
for comparison. 
Observation of the impulse response of the sphere shows 

the following results: 

1. The first peak appears after a transit time across 
one radius of the sphere (which is the time where the inci- 
dent pulse reaches its peak value), and is due to the return 
from the leading edge of the sphere. 

2. The second peak appears at about 5.2 a/C and is due to 
the return from the back side of the sphere, after a transit 


time of: 


Ta 


a, Sted 
C 


Q |p 


which corresponds to twice the free space transit time up to 
the shadow boundary plus the transit time around the back 
of the sphere. This second peak is the creeping wave return. 
3. The response as expected, starts decaying after some 
time equivalent to 5-6 times the transit time across the sphere. 
The ramp response for the same sphere is given in figure 


4.5. The first positive incursion of the response waveform 
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Figure 4.4. Smoothed impulse response of a 


1 meter radius sphere 
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Figure 4.5. Ramp response of a 1 meter 
radius sphere 
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extends to about 5 times the transit time across the radius 
which corresponds to 1.25 round trips across the target. The 
impulse response of a 30 cm diameter sphere was measured in 
the NPS scattering range lab. The waveform obtained is shown 
in figure 4.6. The response appears here with a reversed 
Sign since vertical polarization is used (instead of the 
horizontal polarization considered for the computations). 

The waveform shows very good agreement with the form obtained 
for the 2 meter diameter sphere. 

The second example considered was the cone-sphere whose 
geometry is shown in figure 4.7. The target contour was 
divided into 12 segments while all other parameters remained 
the same as for the sphere. The impulse response of the cone- 
sphere target 1s given in figure 4.8. The waveforms show two 
positive peaks, corresponding to the leading edge and creeping 
wave returns. The separation between them is equivalent to the 
two way free space transit time up to the shadow boundary 
plus the transit time across the rear of the sphere. 

[(2x1.3 +7 x0.7) = 4.7 light-meters] The ramp response 
of the same cone-sphere target 1S given in figure 4.9, and 
the remarks made on the sphere apply here too. 

The third example was the capped cylinder target shown in 
figure 4.10. The contour of the target was divided into 12 
segments, the sampling time was AT = 0.025 m/sec and 
A= 4 oa sec +. The smoothed impulse (gaussian) response 


of this target is given in figure 4.11. The waveform shows 


ae, 





clearly the responses from the leading edge and the rear of 
the cylinder, with the same amplitudes, while from the flat 
horizontal region of the cylinder the response is zero as 
expected from the analytical expressions of the field. 

The ramp response for the capped cylinder is shown in 
figure 4.12. The first positive incursion of the response 


clearly exhibits a shape comparable to the shape of the target. 
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Figure 4.7. Geometry of the cone-sphere 
target (half-target) 
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Figure 4.8. Smoothed impulse response of 
the cone-sphere 
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Figure 4.10. Geometry of the capped cylinder 
target (half-target) 
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Figure 4.11. Smoothed impulse response of a 
capped cylinder (length ~ 0.33 m) 
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Ramp response of the capped 
cylinder (total length ~ 0.33 m) 
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V. INVERSE SCATTERING 


A. PROBLEM DEFINITION 

Given the incident field and the scattered field from an 
unknown target, it is desired in the most general case to 
retrieve the shape, size, orientation, and, perhaps, composi- 
tion of the target. This is known as inverse scattering, or 
target imaging. 

For the present case, the problem is restricted to a 
metallic body of revolution illuminated by a field propagating 
along the axis of symmetry, and the scattered field is measured 
in the backscattered direction, as shown in figure 5.1. 

The solution of this problem is approached here directly 


in the time-domain. 


B. FORMULATION OF THE PROBLEM 

In {6], Kennaugh and Moffatt derive the relationship between 
the physical-optics impulse response of a body and its area 
function A(z), projected in the plane orthogonal to the direc- 
tion of propagation of the field (z-axis). 


The impulse response is given by: 


C A FA 


(Sam) 
as ain Z=Ct/2Z 





P(t) = 


where ry is the range from the observation point to the 
scatterer, as shown in figure 5.1. Integrating twice (5.1) 


with respect to time gives the ramp response as follows: 
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Figure 5.1. Inverse scattering geometry 
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Equations (5.1) and (5.2) were proven to be valid for bodies 
having a monotonic area function with a limiting value at the 
shadow boundary. This relationship between the ramp response 
and A(z) is derived in Appendix F for the body of revolution. 
The total backscattered field that would be measured will 
consist of the physical optics term plus an additional term 
due to the field induced interactions between the currents 
at different parts of the body. The total current induced on 
the body is given by equation (3.1). The first term of the 
R.H.S. of the equation represents the physical-optics current, 
which we will denote by Sy and the second term expressed by 
the integral represents the additional current, which we will 


denote by Jac The total current will then be: 


> > ~ > > > 
J eotay (tr) = Ie ee + Ja (r,t) ( Siusy 


Expressed in terms of (5.3) the total backscattered field for 


a ramp excitation will be: 


S L oh = ca a 
H(t) amr c ™ a-(J_+J,) xa, dA (os } 


_ 3s Ss 
= Bee? + Ha (t) 


where: 


He (t) is the physical optics ramp response, 
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He ae) is the additional response, 


and 


el 


HR 


t) is the total ramp response. 


Using the relationship between the physical-optics response 


and the area function, (5.4) can be rewritten as: 


Heme) = re 50k + Hate) P Ze= Ct/2 as .5) 





which leads to 
Fs s . As 
0 (z) = 2x C (Hp (t) H(t) ) ( 52263) 


Equation (5.6) is the basis for the solution of the inverse 


scattering problem for the body of revolution. 


C. INVERSE SCATTERING SOLUTION PROCEDURE 

Equation (5.6) relates the contour function of the body, 
o0(z), to its ramp response. Solution of that equation leads 
to the shape of the body. The only term known in equation 


(5.6) is the measured value of the total ramp response He (t) 
while Hate is not known. 

Moreover, all the terms in the equation are functions of 
0(z), hence there results a non-linear equation that will be 
solved by an iterative procedure. As a first approximation, 


the unknown term Ha (t) will be neglected and so the first 


estimate of the contour function will be given by: 


o,(z) = [ar oc He (t)) 7/4 ia 


ql 





Once a body is defined its ramp response can be computed since 
we know the incident field. The ramp response computed for 
the first estimate of the body will be denoted by H. (t) and 


iL 
is expressed as: 





2 
mo (2) 
Ss 1 Ss 
meee = FEC + Ho ee) (5.8) 
oO 1 


A new estimate of the shape Pp, (Zz) can be found using (5.6) 


and (5.7), and is given by 


o2(z) = 2x clHS(t) - HR (t)1 + 0? (z) (5.9) 
or equivelently by 
of(z) = 2r Cl2HS(t) - He (t)] (5.10) 
1 


The new estimate p> (2) serves to compute a new ramp response, 
S 
B2 


Following this iterative procedure, the k-th estimate of the 


H_ (z), from which an updated estimate of o0(z) is extracted. 


contour function will be given by: 


2 - 2 Ss _ 
py (2) = Py _y (2) + 2r .C [H, (t) eae (Sr ell | 
Br by: 
Z S “oe S 
py (2) = prod a ales 2 te (+) ) Coney, 


The iterations are continued up to the point where the differ- 


ence between the true measured response and the last computed 


(OY 





response reaches a minimum. Several runs of the program for 
various bodies has shown that, as the iterations go on, the 
mean squared error 1S decreasing, reaches a minimum then starts 
oscillating. In one case observed, the error is diverging 
after reaching the minimum. Some representative results are 
shown in the next section. 

The normalized mean squared error will be defined as 
follows: 


CCS) SBS) 


a2 — ee es keds (5.13) 


S 2 
(Hae) 


2 is the value that must be minimized in order to obtain the 


most accurate representation of the shape of the body. As 
it is defined, the error is equivalent to the difference be- 
tween the successive shapes, as shown by equation (5.11). Once 
the difference He) - He (t) 1S minimized, the difference 
2 2 | —— 
py, (2) - Py _y (2) is also minimized. 
D. NUMERICAL SOLUTION AND EXAMPLES 
For the numerical solution we will use equation (5.12), 
which is more suitable for the algorithm developed. 
We will incorporate the following terminology: 
Hoy (AAT) is the time sample of the measured ramp response 


while He (nAT) is the time sample of the computed ramp response 
k 


for the k-th estimate of the body, and Py. (2) is the k-th 
estimate of the contour function. In discrete numerical 


form, equation (5.12) becomes: 


ig 





Kaa 


2 —— s —_ 
eee cee Hy, (AL) H, (nAT) | (5. 45 
2=0 i 
where 
Zz, = 4 CT 72 
aC mEOle. = Ui, He (nAT) = 0. The mean squared error is: 
e) 
e 2 
) (Hay (nT) Sis (nAT) ] 
2 n=1 k-1 
EL a = GSS } 
2 
) Hee, (282) 
n= 


where N is the total number of time samples observed and 
computed. 

The ramp response of the body defined by Py, _3 (2,) will 
be computed using the integral equation solution to the for- 
ward scattering problem, as already described in Chapter IV. 

The Fortran program, given in Appendix E, has been extended 
with a subroutine named "SHAPE", in order to implement the 
inverse scattering algorithm. In this subroutine the shape 
of the body is computed iteratively, using (5.14). 

The space samples given by (5.14) are defined at equal 
increments on the z-axis, while the program developed for the 
computation of the ramp response has been proven to work best 
when the contour was divided into equal length segments. Con- 
sequently, one of the tasks of subroutine "shape" is to find 


the space-samples at the points delimiting segments of 
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approximately equal length. This is done using a "cut and 
try" method which converges to the desired length of segment. 

The iterative approach for solving the inversion equation 
(5.14) was tested in a simulative way for the three types of 
bodies already considered for the forward scattering solution 
in the previous chapter. 

The ramp response of the body is first computed and is 
used to simulate the response thatwould have been measured 
(without noise or error). Then the shape of the body is re- 
trieved from that ramp response. 

The shape of the body is retrieved from the first positive 
incursion of the ramp response, which as it will be seen 
extends beyond the length of the target for the shapes 
considered. 

The first example was the 1 meter radius sphere. Figure 
5.2 shows the l-st, 3-rd, 6-th estimates of the contour func- 
tion of the sphere. The 6-th estimate gives the least error 
and is shown superimposed on the original shape in figure 5.3. 

| The second example was the cone-sphere, used in two con- 
figurations. The first one when the contour of the body was 
Originally cut to 12 segments. Figure 5.4 shows the successive 
estimates (1,3,6) of the contour, and figure 5.5 shows the 
Minimum error estimate (estimate 6). The second version, 
when the body was cut to 1l segments, gave a much more accurate 


result, as shown in figures 5.6 and 5./7. 
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The last example was the capped-cylinder. Figure 5.8 
shows the successive estimates of the shape (estimates 1,3,4), 
and figure 5.9 shows the 4-th estimate which gives the minimum 
error. 

For the case of the capped cylinder, after reaching the 
Minimum error at the 4-th iteration, the imaged shaped of 
the body started to diverge considerably, and appeared more 
shortened than the actual shape. 

Common to all examples, the error appears mostly at the 
shadowed region of the body, but in general the shape is 
retrieved with quite a good fidelity, which proves that the 
iterative direct time domain solution is working, subject to 


certain caveats which will be discussed in the conclusions. 


EB. INFLUENCE OF NOISE 

In the examples considered in the previous section the 
only noise appearing in the ramp response is the numerical 
noise generated by sampling and computational round-off. In 
this section the measurement noise will be tested with a ramp 
response to which noise is added. In this simulation, each 
noise sample is obtained by the sum of a large number of zero- 
mean, uniformly distributed random numbers, which, by use of 
the central limit theorem, approximate the statistics of 
gaussian noise. The noise samples are normalized in order to 
remain within a certain voltage signal to noise ratio, and 
then are added to the computed ramp response. 

The body considered for the noise test was the sphere, 


and the signal to noise ratios tested were: 
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V 
20 log = = 10, 20, 30 and 40 4B. 


noise 


The sphere ramp responses with the noise added are shown in 
figures 5.10, 5.11, 5.12 and 5.13 for the successive signal 
to noise ratios. 

For the 10 dB S/N case, the appearance of spikes and dis- 
continuities on the ramp response causes the shape of the 
body to diverge considerably since the first iteration, and 
the shape cannot be recovered. 

For the 20 dB S/N case, the computed shape starts to con- 
verge, with the minimum error shape obtained at the third 
iteration. Then the imaged shape proceeds to diverge and to 
shorten. 

Figure 5.14 shows the successive estimates of the shape 
up to the 3-rd, and figure 5.15 shows the sixth estimate, which 
is clearly different of the actual shape. 

For the 30 dB S/N case, the minimum error shape was ob- 
tained at the 5-th iteration, as shown by figure 5.16. Then 
again, the computed shape diverges as shown in figure 5.17. 
Common to all cases of divergence, the wrong estimates start 
as soon as the last computed shape exhibits a strong discon- 
tinuity in its contour. A very important conclusion from the 
noise analysis is that before starting the inverse scattering 
procedure the ramp response used must be smoothed so that no 
spikes or sudden discontinuities appear on it. 

We should normally expect that the true measured response 


will have a 20 dB S/N ratio, so that the problem of divergence 
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Figure 5.10. Ramp response of the sphere with 
noise added (S/N = 10 dB) 
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Figure 5.11. Ramp response of the sphere with 
noise added (S/N = 20 GB) 
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Figure 5.12. Ramp response of the sphere with 
noise added (S/N = 30 dB) 
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Figure 5.13. Ramp response of the sphere with 
noise added (S/N = 40 dB) 
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must be investigated. In any case, before diverging the 
shape that gives the minimum error is quite a good approxi- 
mation of the true shape. 

An additional test with 40 dB S/N was made, and no signi- 
ficant deviation, from the case without noise added, was 
observed. Figure 5.18 shows the first and 6-th estimates of 
the shape. The minimum error shape was the 6-th estimate. 

Up to the 12-th iteration the shape of the body remained 


stable, with minor variations. 
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VI. SUMMARY 


A. SUMMARY OF RESULTS 

The inverse scattering method for radar target imaging 
was the main topic addressed in this work. The inverse 
scattering algorithm developed uses the time domain solution 
of the integral equation. The class of targets considered 
included the conducting bodies of revolution whose important 
symmetry property reduces considerably the computation time 
and the computer storage needed for the solution of the in- 
verse scattering problem. 

The use of the ramp response for retrieving the shape of 
the target proved to be a very efficient tool, and the exam- 
ples considered gave generally good results, with initial 
convergence of the imaged shape. For the case of the long 
cylinder and the sphere with low signal to noise ratios there 
appeared an eventual divergence of the imaged shape after 
many iterations. This is thought to be due to accumulated 
numerical roundoff errors in repeatedly solving the integral 
equations. In general it has been shown that the use of 
the ramp response by itself, without any further computation, 
gives a close approximation to the shape of the body but up 
to its shadow boundary which corresponds to the maximum of the 
area function A(z) along the line of sight. Beyond that boundary 
the body is, at the first observation, elongated and the 


iterative algorithm is needed to compensate for the errors 
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induced by the wake or "creeping wave" which alters the 
target response. 

As the iterative computation process goes on, some numeri- 
cal instability that is mainly due to a discontinuity in the 
computed shape might arise. However, before diverging, the 
Shape of the body always reaches an estimate which has a 
minimum error relative to the true shape, and is a good 
approximation of that shape, so that the minimum error cri- 
teria set for stopping the iterative process appear to be an 
attractive way to overcome the instability phenomena. As a 
continuation of this work it is suggested to carry out a 
detailed analysis of the error and numerical instabilities 


that have been observed while testing the iterative algorithm. 


B. RECOMMENDATIONS 

The algorithm developed for target imaging by transient 
inverse scattering, has been tested numerically by simulating 
additive physical noise. It must now be evaluated in a real 
laboratory environment where the true measured response is 
used. The transient scattering range at NPS should be used 
to test and finalize the algorithm for the class of targets 
addressed in this thesis. 

The first step to be done is to obtain the measured ramp 
response of the target. The range uses an impulse generator 
which exhibits a gaussian pulse shape characteristic. The 
true ramp response can be obtained by using either a time 


domain convolution or a frequency domain computation and I.F.T. 


oF, 





The true ramp response can be obtained, in principle, from 
the gaussian response of the target using the following 


relationship: 


H (w) 
_ poe ee 
Hptw) = Glo) 3 


where: 


Hp (w) 1s the Fourier transform of the true ramp 
response of the target, 


H..(w) 1s the gaussian impulse response of the 
IG 
Cart gcicy 


G(w) is the Fourier transform of the transmitted 
gaussian impulse, 


and the ramp response will be: 
Ba(t)) = F “AH (w) | 
R R 


Very special care must be given to the choice of space 
sampling and time sampling intervals. The time sampling is 
dictated by the sampling rate of the oscilloscope, and the 
Space sampling must be done at a high enough rate such that in 
every case the distance R between the closest spatial samples 


1s related to the time sampling T by: 
Ry, Ce AT 


The next step is to translate the FORTRAN program given in 


SDOSt el i eon ciceCcOMPUteCrE Ob the labw There isea need 
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to reduce to a minimum the computation time since the 
laboratory microcomputer is very slow compared to a large 
Mainframe. The use of some additional symmetries of the 
body can certainly help in this task. 

Finally, it 1s necessary to Minimize the noise in measure- 


ments, and smooth the response to be used for shape computation. 


C. AREAS FOR FURTHER STUDY 

The ultimate objective of this research program is to 
aid in determining the feasibility of developing future radar 
systems having the innate ability to classify targets using 
their transient scattering response. 

In the most general case the target will be of arbitrary 
Shape and the use of the iterative method described in this 
thesis would be highly impractical and prohibitive for com- 
plex targets, due to two principal factors. The first is the 
time needed to compute the response of a complex body, and 
the second is that the area function obtained from the ramp 
response 1s not sufficient to determine in what way the 
response 1s to be computed. 

The best technique to deal with such targets would involve 
a way that does not require the computation of the integral 
equation, and instead will use only the measured impulse 
response (or ramp response) to discriminate between targets. 
The first solution involves measuring the ramp response at 
multiple look angles. The response waveform will be used up 


to the shadow boundary to find the area function along each 


oie, 


line of sight. The target image can be approximately recon- 
structed from the projections of the area functions. Das 

and Boerner [18], suggested a very interesting approach for 
target shape estimation using these projections. This method 
seems attractive and is felt that it could lead to a good 
solution for target imaging. Its main disadvantage is the 
need for multiple look angles, requiring either interfaced 
multiple radars, or long-time dwelling with one radar. 

An alternative approach is to classify targets via the 
complex poles and residues of their time-domain signatures. 
Once the poles and residues are extracted they could be com- 
pared with stored library data corresponding to Known targets, 
and thus would provide a classification of the target under 
test. This catalogue approach is also being pursued using 
the NPS transient range, as well as through other independent 
efforts, and may prove to be the most practically implementable 
technique available. 

Another area which could form the subject of future re- 
search includes the time correlation of the shape of the 


target with its motion. 


100 





APPENDIX A 


DERIVATION OF THE SCALAR INTEGRAL EQUATIONS 


In Chapter III, equations (3.3) and (3.4) involve some vector 


multiplications. The first appearing in the integrals is: 


w~ 


a xd'xR = (a xa' xR)-Jd' + (a xa! xR)g! (Al) 
n ig) S S n oO) vu) 
We note the following identities: 
a xa'xR = a'(a -R) - R(a.-a!) (A2) 
n S s on n s 
A a > = oar a = — “A cae 
an ACY xR avon R) BME oat.) (A3) 


F > “Aw “~ p “~ “aw “~ . 
Expressing R, a' and a! in terms of the fac7a ray) coordinates 


S p 


system gives. 


R = a, (Rea) “- a, (R-a.) 
a 4 = ay (agray) + ae \44° 5? (A4) 
a = a, (ac-a,) + a, (al-a,) 


For the body of revolutions we note the following relationships: 


Pa. an 
a'-a = cos 8 
d 9 
alea = - sin a-sin 8 
d Ss 
ras wn A 
a'ea = - cos a-sin 8 
dg n 
aca, = sin a sin a'tcos @ + cos a cos a' 


Oa: 





acca == Siiewmeasin 6 


eosmasSim @'cos 6 — Sin @ cos q' 


0) 
ll 


where 


> “aw “a “aw ; : F 
The vectors R, a_, a_, a, are given in carteSian coordinates 


Sy 2 
R = a,(p cos¢-p' cos $') + a (o sind -p' sing") 
7s too 
a, (2 Zz) 

a = a (cosa-cos od) + a (cosa-sin $) - a sina 
n < 4 Z 

a = a (Sina-cos $) + a (Sina-sin ¢) + a cos a 
S x y Z. 

a = ~a (Sin a) + a _-cos 

p x : Y \ 


Carrying out the dot products gives: 


Rea = cos a(p -op'cos £8) + (z' -z)Sin a 

Rea. = sin a(p -p'cos 8) - (z' -z)cos a (A5) 
> «A 2 

ae = - p' sin 8 


Using the relationships given by the sets of equations (A4) 


and (A5), and simplifying, equations (A2) and (A3) become: 
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> A 
4 . rR a, isin B{(z'-z)sinasina' + pcosasina' 


= o'cosa'sina]} +a_{(p-p'cos8)cosa' 


+ (z'-z)cos8sina'} (A6) 
a al xR = a, ((2'-z) sinacos B+ cosa (opcosf-o') | 
+ ely samieliz )—Z) (A7) 


Defining: 


Ve 7 — BZ 
Fy = (z'-z)cos ®8sin a' + cos a'(o-o'cos 8) 

U = (z2' -z)sinasin o' +ocosasin a' -o' cosa'sin a 
FE. = cos a(pcos 8-9') + sin a cos 8(z' -Z) 


produces from (A6) and (A7) 


a x4' xR = 4 U sing + a -P (A8) 
n S d S ae 
a, * ay xR = a, m= a. V-sin 8 (AQ) 


Note that V, Poo FE. are even functions of 8 and U is inde- 
pendent of 8. Next we have to carry the second vector multi- 
plication: 


a xH = -a_H’ sin 6 - a, H’ cos 6 sino (Al0) 


Y 
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me 
| : : : : 
VicrChime ls ag LVCim Lae cabtesian Coordinates by: 


After substituting equations (A8), (A9), and (A100) into 
equations (3.3) and (3.4) we obtain the scalar expressions 
for the components of the current in the S and 6 directions 


as follows: 


op eel: iy i ee an . 
ia sciesin ¢ + d-slete aa Lepas J, Vsin B]dA' (All) 
ee 1. I gel 
J, = 72H cos ¢sina +> | =glle oe Fe. 
AR 
° U j ' ' 
[J2 u sin B + TE oidA (Al2) 
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APPENDIX B 


DEMONSTRATION OF THE SINUSOIDAL VARIATION 
OF THE CURRENTS WITH 96 


Maxwell's first equation: 
> 
* a 
VxH = ce E+ (Bde) 


In the case of interest, due to symmetry, the field H is 


independent of » and is given by 


=} A 
H = 2. He ,2,t) 


or in cylindrical coordinates 


> A ; A 
H = A cos 0 as = Hy Sin 06 an 
dH 
a i 9 ; ia O a 
VxH = LS 5g (0 HJ isin ¢ a, + 3777608 ba, 
oH 
A, Sia: (B2) 
a. Sin o 5p 


The currents on the body expressed in cylindrical coordinates 


are: 

-> “Aw ‘ “~ “A 

J = a SinaJ._+a,JI, +a cosa J (B3) 

Pp S d © 72 S 

The electric field: 

Ee=s a El = 3 E +4 E (B4) 

= a _ a Sin 9 4 a cos és 

Substituting equations (B2), (B3) and (B4) in Maxwell's 


equation (Bl) and equating the terms in the same direction gives: 
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, s As - | 
Sin a J. = a yq (PH) Eo]sin d 


The term inside the brackets is independent of 9» so that we 


can write 


Siied B=8) (0,Z,b)sin. ¢ (B5) 
S © 
iL 
Similarly 
OH A 
ao = aA. = E,)cos 0) 
Je = shes ee cos 0 (B6) 


and finally 


JH 


O 
JI cos @ = - sina =—-—— 
S dp 


J  (9,2,t) sin ¢ (B7) 
we 


Equations (B5), (B6), (B7) show clearly that the currents 


J and J, vary sinusoidally with the angle ¢ so that we can 


Ss ® 


generalize: 


J(0,2,t,9) = J(9,2,t) sin ) 


J, (0 ,2,,9) — Ty Wo oie: cos 
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APPENDIX C 


SIMPLIFICATION OF THE SCALAR INTEGRAL EQUATIONS 


In Chapter III the scalar equations for the currents 
components were given by equations (3.4) and (3.5). In these 


equations we had: 


dA' = op' ds! dé 
and 

JS = Js Sin 0 

Jo = mM cos @ 


Substituting in equation (3.4) gives 


. 27 
, i ne se ieee. 1 
Js Sin Oe ee 2Hesin @ + a J p'ds i ale 37 +R! 
0 0 R 
-[Je sin vB os Sin 8 Vcoso'jdp (C1) 


As notedin Appendix A, the expressions V, FS and F. appearing 
in the integral are even and periodic functions of 8. We 
recall that 

LS (07 +p'* - 2ppcos 8 ae See 


and 


T = t - R/C 
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so that R, Je(P',2,T) and J'(p',z',t) are also even and 
periodic functions of 8. Finally Uis independent of 8. 
We will concentrate on the integral over 8 in order to 


reduce it to its simplest expression. Since 


we can write 


Sin 9' Sin @ cos 8 + sin 8 cos 6 


cos ¢' = cos @ cos 8 - sin » sin 86 


and substitute this into the integral over 8 in 

equation (Cl). The integral over 8 for a complete period, of 
all terms that contain sin 8 will be zero since those terms 
are periodic and odd functions of 8. All the other terms will 
produce a non-zero integral and will be retained. 


After this simplification, equation (Cl) will become: 


; Ss 
] = = + - : IL 1 , 
J. Sino? = 2H” sin > + sin 6 oT ) o'ds 
zal 
e 1 0 J. 1 i] : 2 
y) = +t +=] [JiFS cos 8 ra Gabe RjJde (C2 ) 


which after repressing the common term sin 9, leads to the 


final equation: 


; S 27 
1 i 1, ak @ 1 
J = -2H + = f p'ds' i sla a +S] 
Ss 27 0 0 R Cote RB 
ee 
e 1 1 
[JSP cos 8g +VI,sin RJdé. (C3) 
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Following the same reasoning for equation (3.5), and repressing 


the common term cos ¢ leads to the following final equation: 


f J (* _" iE 
a) = -2H sina+t+ += ods. —[ 
db aT 0 0 R 


,) ie 
37 +R! 


Ole 


“(gt U Fin-g et ot 


iF 2 cos 8]d8 (C4) 


In both equations (C4) and (C3) the currents are: 


J sed 7 ee a 


ae ae 


t being the retardated time t - R/C, and the incident field 
1s: 


H> = 4H(z,t) 
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SYSIN DOD DATA 


APPENDIX D 


ORE AE SC A eC RE a a RC aK a CK 2 KC BK a 2 2K RE a eK ak eo a a a OR a a Ae a oe 
a * 
é TRANSIENT RESPONSE OF A STRAIGHT THIN WIRE a 


BE EAE He AEH AK He CK HE EA A HG EE AS RCT he eA ee AE 2 He He He ee ae he He he Ke Ae OK 


PURPOSES: 
Foes Gap mg 

THIS PROGRAM COMPUTES THE TRANSIENT TIME DOMAIN 
RESPONSE GF A STRAIGHT THIN WIRE. 

THE PROGRAM COMPUTES THE CURRENTS INDUCED ON THE 
WIRE BY AN INCIDENT IMPULSE FIELO,AND THE SCATTERED 
FIELD, FOR ANY ANGLES QF INCIDENCE AND SCATTERING. 

USAGES: 
Me ek He 2k 4K 

THE PROGRAM CAN BE USED AFTER PROPER DEFINITION 
OF THE INITIAL PARAMETERS QF THE WIRE,AND DEFINITION 
OF THE DESIRED ANGLES OF INCIDENCE AND REFRACTION. 

THE WIRE [S ALIGNED WITH THE Z-AXIS STARTING AT Z=0 
THE ANGLES ARE MEASURED WITH RESPECT TO THE Z-AXIS 

PARAMETER DEFINITION 

He Beste ake Heake ake oe whe te he we Heke He Xe Se ak ae x 

BETA(400740) :ARRAY DEFINING THE CURRENTS FOR UP 

T2 40 SEGMENTSAND 400 TIME INCREMENTS.) 

GAMMA{ 400,40) ARRAY DEFINING THECHARGES FOR UP 

TO 40 SEGMENTS AND 400 TIME INCREMENTS, 
NS: NUMBER OF SEGMENTS IN WHICH THE WIRE IS DIVIDED 


SLs TOTAL LENGTH OF THE WIRE IN METERS. 
Az RADIUS OF THE WIRE IN METERS. 
ITS: NUMBER JF _ TIME INCREMENTS. 
pang canoes A MULTIPLE OF 80 FOR PLOTTING 
TI: ANGLE OF INCIDENCE IN DEGREES. 
TS: ANGLE OF REFRACTION IN DEGREES. 
DZ: LENGTH OF A SEGMENT (L/NS) 
DTs TIME ENCREMENT IN SECONDS (O0T=D2Z/C) 


AAs GAUSSIAN IMPULSE PARAMETER.IN SEC-1l. 
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1 GAMMA (400 140) ,ES( 400) 


ix 
3 
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MAIN PROGRAM 
HK RH KK KH KE KK 


MK Ye HE Ree HC He ee eke HE HK HK 
* 
# 
$e 


OVUVOOOUYUOU 


CONSTANTS INITIALIZATION 


WOOHOO 


PARAMETERS INITIALIZATION 
COMPUTATION PARAMETERS 


WIRE PARAMETERS 


QUO 


E OF INCIDENCE TI 
oO, 999 
DEGREES TG RADIANS. 
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GO 


Od T1,TS 
WRITE HEADING 


GAUSSIAN IMPULSE PARAMETERS 


TRANSLATION FROM 


CONTINUE 
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340 


330 


OOM 


WOOC 


o) Eid 


333 
33% 


350 
300 


MAONANMONONDA OMOOM 


OOO 


OM 


0) GO TO 330 


Nr" 
ene @ mM 


Ooorc 
mA TUTTO 


Saree £5371) 


meOOO NIMS 
mM 


namOom 

AUS UpPAK KH 
mt TT ed OL PO RIND mt 
emmeaQoite & 


arar it tta2 


LE.90 
FOR ANGLES OF INCIDENCE MORE THAN 90 DEGREES 


L2=L-NS 
GO TO 332 


FOR ANGLES OF INCIDENCE LESS THAN 90 DEGREES 


CONTINUE 

ee 

CONTINUE 

iF (J.GT.50 ) GJ TO 333 

El= EXP (-AA¥*2) *( SOT #(1-L2* AT )-TMAX) ®*2) *SIN(TIR) *3 
GO TO 334 

CONTINUE 

EI=0.0 

CONTINUE 

BETACJ»L) = EL -DELAZ /OT-DELPHI /DZ}/(CL#FO/0T) 
Li=t-1 

[F(C.GT.1) GO TO 350 

GAMMACJSyLI= GAMMACJ1,L I-BETA(JyL I/C 

GO TO 300 

CONTINUE ; 

GAMMACJyL)= GAMMA(J1yL I-BETACSsL 9/0 #BETACUyLLI/C 
CONTINUE 

GAMMA(JyNS)= GAMMA JL» NSI#BETA(JSeNSLI/C 

KK=NS/2 


WRITE THE CURRENT INOUCED AT THE CENTER OF THE WIRE 


WRITE(5e¢9690) JeBETAC Js KK) 
FORMAT(1OX,14,21X,E12.5) 
CONTINUE 


END OF CURRENTS COMPUTATION 


START CF SCATTERED FIELD COMPUTATION 


0G a0 LP=1,1TS 
ES(LP)=0. 
DO 500 K=1,NS1 

If (TS.GT.90-0) GO TO 501 
IF (TS.EQ. 90.0) GO TO 503 


FOR REFRACTION ANGLE LESS THAN 90 DEGREES 
NQ=UP=-1-(KP-0.5 


NQP=LP—1- (KP+0 
GO TO 502 


»*AR 
5) *AR 
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FOR REFRACTION ANGLE MORE THAN 90 DEGREES 
FOR REFRACTION ANGLE EQ. 90 DEGREES. 


GO TO 502 
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(LP) +DZ*(BETACNQ,KI-BETA(NQL,K)) 


(LP) +DZ*2* BETAI(NQsK) 
S(LP)+4+DZ*¥2* BETACNQP,K) 


eNO 


PLOT OF CURRENT AY CENTER Of THE WIRE 
PLOT EVERY 4M POINT. 


END OF FIELD COMPUTATION 
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PLOT OF SCATTERED FIELD 
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END OF PROGRAM GO BACK TO NEXT ANGLE IF ANY 
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EXAMPLE OF DATA CARDS 
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END OF PROGRAM FOR INVER=0 AND INVER=1 
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IF {INVER .£Q-1) GOTO 19 
CONTINUE 
STARTING POINT FOR INVER=2 
iINC=2 
READ THE GIVEN RAMP RESPONSE 
DO 12 N=1,ITS 
READ (559520) HSEM (N} 
FORMAT (1F10.5) 
HSE (NJ =HSEMON) 
CONTINUE 
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CONTINUE 
STORE THE COMPUTED RAMP RESPONSE 
D0 910 N#1,ITS 
HSE M(NI=HSE IN) 
CONTINUE 
CONTINUE 
STARTING POINT FOR INVERSE SCATTERING COMPUTATION 
ITER=1 


COMPUTATION OF THE FIRST APPROXIMATION OF THE SHAPE 


CALL SHAPE (OSeITER CP ySO0rIT) 
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CONTINUE 
COMPUTE THE RESPONSE GF THE APPROXIMATE BODY. 


CALL FIELD (DS,CP,yDT,INC) 
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79666) EPS 
{ 1 EAS 0 EE oS MEAN SQ. ERROR® ,1F1O.5) 
eGT.1) GO TO 961 


END THE PROGRAM WHEN THE ERROR IS MINIMUM 


IFCEPS~GT.eEPS1) GO TO 999 
EPSISEPS 


CONTINUE WHEN THE ERROR IS NOT YET MINIMIZED 
ITER=ITER+1 
EXTRACTION OF THE PHYSICAL OPTICS RESPONSE 
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E PHYSICAL OPTICS RESPONSE TO COMPUTE A NEW SHAPE 
CALL SHAPE (OS;,ITER,CP ,DSO,IT) 


RETURN TQ COMPUTE THE RESPONSE TO THE NEW BODY 


Go TO 959 

CONTINUE 
END OF PROGRAM WHEN MINIMUM ERROR ACQUIRED 
WRITE FINAL ERROR 

RETEMEG, 96259 LTER+=PS1 

ORMAT (//;20X,®EINAL SHAPE AFTER", 14,"ITERATIONS® /// 
+20X,"NORMALIZED MEAN SQUARED ERROR IS?,1F10.5) 
CONTINUE 
END OF PROGRAM 

STOP 

END 
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GAUSSIAN IMPULSE PARAMETERS 
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TMAX=DT*5.0 
AZ =AA*AA 
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WHEN INCIDENT FIELD IS A RAMP 
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CONTINUE 
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APPENDIX F 


RELATIONSHIP BETWEEN RAMP RESPONSE AND SHAPE 


The validity of (5.2) will be demonstrated analytically 
here in the case of a metallic body of revolution. We first 


express the axially directed incident ramp field as: 


| eo 27 e tee 
2: — 
Hp (2,t) = er) 
0 (ace 2/6 


The physical optics current induced on the body is given by: 


— 


SS “a >1,> 

J(r,t) = 2a, x HR (r,t) Oe) 
The backscattered far-zone field is given by equation (3.8) 
and its scalar equivalent for the body of revolution is 
expressed by (3.16). Substitution of the current in (3.16) 


gives the physical optics ramp response as follows: 


1 
4ir C 
fe) 





Ss : 
J S—[-4H>(z',t)sin a']p’ds’ (F.3) 


The time derivative of the incident field expressed by (F.1) 


is unity: 


i. 
oH, 
Ot at 


1B te <=27C (F.4) 


Substitution of (F.4) into (F.3) gives: 
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: : (* sin oO! 6! ds’ (F.5) 


From the geometry of the body of revolution we note that: 
do’ = sin a’ ds’ 


Substitution in (F.5) gives the following expression: 


i ae a2 


S -_ ' ' = — 
Bye) SM Sf uci? ae mec & (2) (F.6) 


Since the area function is given by: 
ANZ.) 5 aT 9*(z) 


the physical optics ramp response (F.6) will be related to 


A(z) by 


t) = A(z) Ca) 
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